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The properties of the strongly interacting edge states of two dimensional topological insulators 
in the presence of two particle backscattering are investigated. We find an anomalous behavior of 
the density-density correlation functions, which show oscillations that are neither of Friedel nor of 
Wigner type: they instead represent a Wigner crystal of fermions of fractional charge e/2, with e 
the electron charge. By studying the Fermi operator, we show that the state characterized by such 
fractional oscillations still bears the signatures of spin momentum locking. Finally, we compare 
the spin-spin correlation functions and the density-density correlation functions to argue that the 
fractional Wigner crystal is characterized by a non trivial spin texture. 
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The helical Lnttinger liquid[l] (hLL) is the state of 
electronic matter that describes the interacting helical 
edges of the recently predicted [2-4] and experimentally 
realized[5, 6] two dimensional topological insnlators[7- 
9] (2DTI). Particnlar emphasis has been devoted to the 
investigation of the transport properties of 2DTI: topo¬ 
logical protection of the edge states facilitates the ob¬ 
servation of condnctance quantization[5, 6, 10] in short 
samples. Long edges on the other hand are characterized 
by a reduced conductance. Even though a comprehensive 
theoretical understanding of the scattering sources caus¬ 
ing the reduction of the conductance of the edge is still 
lacking, the role of electron-phonon interactions]!!], of 
magnetic[!2] and nonmagnetic impurities[!3, !4], in the 
presence of random Rashba disorder[!5, !6], of break¬ 
ing of axial symmetry in combination with electron- 
electron interactions and impurities]!?, !8], of tunnel¬ 
ing among the edges and charge puddles in the bulk 
of the 2DTI]!9], and of the coupling between opposite 
edges]20, 2!] has been theoretically elucidated. The 
mathematical tool allowing for most of such calculations 
is bosonization]22, 23], a procedure that enables us to 
recast the Hamiltonian of the interacting electrons on 
the edges into a Hamiltonian of free bosonic excitations, 
representing charge density waves, and to express the 
Fermi operator in terms of the creation and annihilation 
bosonic operators. The physical meaning of the bosoniza- 
tion technique can be understood within the framework 
of Luttinger liquid theory ]24], that is the one dimensional 
counterpart of the Fermi liquid theory for one dimen¬ 
sional gapless systems. More precisely, the exactly solv¬ 
able Luttinger model[25-27], a strictly linear theory of 
interacting one dimensional fermions with infinite band¬ 
width in the single particle dispersion is usually em¬ 
ployed. The validity of the Luttinger model as a basis 
for the description of interacting electrons has a num¬ 
ber of experimental demonstrations, ranging from spin 
charge separation]28], to charge fractionalization]29, 30], 


and to anomalous tunneling[3!, 32]. On the other hand, 
the Luttinger model alone fails in predicting a reason¬ 
able behavior of local observables[33-36], such as the 
electron density and the density-density correlation func¬ 
tions when electron-electron interactions are strong. In 
particular, the Luttinger model is not able to capture the 
transition between a weakly correlated state dominated 
by Friedel oscillations of the density [37, 38], and the 
strongly correlated one dimenional Wigner crystal [39- 
46]. To overcome this problems one has to consider a 
richer theory: the Luttinger liquid, the universal model 
describing low energy properties of gapless one dimen¬ 
sional systems [22-24]. The construction by Haldane [24] 
clearly shows that the Fermi operator ifsix), where s = ± 
is the spin projection, of a generic one dimensional elec¬ 
tron system can be significantly different with respect 
to the one of the Luttinger model[24, 36]: the stan¬ 
dard relation if.ix) - with 

kp the Fermi momentum and 4i{x) and 9{x) the usual 
bosonic fields, is replaced by the more general expression 
ips ~ with the model 

dependent coefficients 

We aim at understanding the properties of the strongly 
interacting hLL, and hence we have to build the appro¬ 
priate Luttinger liquid theory. This formulation presents 
difficulties since some of the usual paradigms break down: 

{i) spin momentum locking breaks the symmetry, usu- 

(s') fs) 

ally holding for standard Luttinger liquids, c_^ = Cp, 
since it would imply no preferred chirality for a given 
spin projection, 

(m) time reversal symmetry protects from one particle 
backscattering off nonmagnetic impurities, which implies 
that the usual Friedel oscillations of the density are for¬ 
bidden: in fact, if they were present a capacitive coupling 
between the electrons of the edge and the impurity poten¬ 
tial would lead precisely to one particle backscattering. 

In this Letter, we develop a Luttinger liquid picture 
of the strongly interacting quantum spin Hall system in 
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the presence of two-particle backscattering extending 
over all the helical edge. We find a state characterized 
by charge oscillations. These oscillations are profoundly 
different from usual Friedel or Wigner ones in view of 
their different wavelength: they are characterized by a 
wavelength that is half of the wavelength of the usual 
Wigner crystal, suggesting the formation of a correlated 
state of fermions with charge e/2, with e the electron 
charge. Moreover we show, by studying the series 
expansion of the Fermi operator, that peculiar features, 
inherited by spin momentum locking, still survive in 
the strongly interacting regime characterized by density 
oscillations. Finally we address the spin-spin correlations 
and demonstrate that they have the wavelength of the 
usual Wigner oscillations, suggesting a complex spin 
pattern. 


The classical solution i/o of the equation of motion, upon 
which a low energy theory for fluctuations can be ob¬ 
tained with standard techniques, becomes [50] 


(j)o{x,X) 



m{x + X) j 
k 


( 4 ) 


with am [y, k] the Jacobi amplitude function with ellip¬ 
tic index k. The quantity X € [0,L] is a free param¬ 
eter whose presence is due to translational invariance. 
Physically it represents the ’center of mass’ of the soli- 
ton solution. The index k is fixed by the the number Q 
of electrons on the helical edge, as measured from the 
Dirac point: since the density po(x), corresponding to 
the classical solution, is given by 

po(x,X) = ^d^(/>o(x,X), (5) 


The Hamiltonian of the system in its bosonized form 
can be written as 


H = 


dx ■ 


n 

2 




- + h2p \ ( 1 ) 


where the first contribution is the usual Luttinger liquid 
Hamiltonian, the second is the chemical potential term, 
and h 2 p, given by 


h 2 p = cos(/3(/)), 

z(7raj^ 


( 2 ) 


is the two particle backscattering. More specifically (p is 
the Luttinger liquid bosonic field, ]/[ its conjugated field, 
p is the chemical potential, /3 is related to the Luttinger 
parameter Kl hy j3 = Jy'FiFi, g 2 p measures the strength 
of the two-particle backscattering and a is the cutoff of 
the Luttinger theory. Periodic boundary conditions, with 
period L, are imposed. Two-particle backscattering can 
arise in the presence of anisotropic spin interactions[l] 
or in generic helical liquids [47]; it is allowed by time re¬ 
versal symmetry and its effect is more pronounced when 
the Fermi level is close to the Dirac point. Our interest 
in such an interaction term is due to its formal analogy 
with the umklapp term occurring in usual Luttinger liq¬ 
uids, which is known to lead to Wigner oscillations in the 
density [22, 33], and finally to the formation of Wigner 
molecules[48] in finite systems. 

We are interested in analyzing the helical counterpart 
of the Wigner crystal, and hence we focus on the regime 
of strong interaction /3 1. An effective strategy to 

deal with such a regime in the case of a regular spinful 
Luttinger liquid has been developed in Ref.49. The same 
strategy can also be applied to the case of the hLL. The 
first step, after defining for notational convenience the 
mass m = v^52p/3^/(27r^a^) of the so-called breather, is to 
recast the problem into the Lagrangian formalism. The 
Lagrangian density C reads 

{dtpf {d^4>f io\ 

^ ^-^- yy cos{l3(p) + -^dx4>- (3) 


one has 


mL 

'Yq 


kK{k), 


( 6 ) 


with K{k) the complete elliptic integral of the first kind. 
In fact, by integrating Eq.(5), we obtained (j){L — X, X) — 
(j){—X,X) = AirQ/fd; the same difference can be com¬ 
puted using Eq.(4). By matching the two expressions 
one obtains the result in Eq.(6). 

The number of particles Q is instead controlled by the 
chemical potential, via the minimization of the total en¬ 
ergy. The condition expressed in Eq.(6), which is differ¬ 
ent from the one occurring in usual Luttinger liquids[49], 
has important implications on physical observables. Al¬ 
though a deeper analysis, which will be discussed below, 
requires the inclusion of quantum fluctuations, a flavour 
of the physical consequences of two-particle backscatter¬ 
ing can be gained already by considering the properties of 
the classical solution. We can consider the electron den¬ 
sity po{x,X) once a particular choice of X, say A = 0, 
is done. This quantity describes the local electron den¬ 
sity if a local perturbation able to pin it is introduced. 
It is important to note, however, that when the aver¬ 
age (integration) over X is carried out, as required for 
the clean system described in Eq.(l), the average local 
density po = QjL is constant due to translational invari¬ 
ance. Alternatively we can consider the density-density 
correlation function ^{x) in the classical regime, after av¬ 
eraging (integrating) over the variable X. The average 
over X is physically needed since the center of mass of 
the classical solution can be anywhere in the ring with 
equal probability. Explicitly one has 


2QKik) 

Pq{x,0) = --— dn 

ttL 


4:QK{k)x 


,k 


^o{x) = 


1 


dXpoix,X)po{0,X), 


( 7 ) 


( 8 ) 


where dn [a;, k] is the Jacobi dn function with elliptic in¬ 
dex k. As shown in Fig.l, the number of peaks of po{x, 0) 
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and ^(a;) is 2Q, instead of Q, as it is in the usual Luttinger 
liquid, and as would be expected if the system was in the 
Wigner molecule regime. It will be soon shown that this 
feature is also present when quantum fluctuations are 
included, and is hence a characteristic of the helical liq¬ 
uid: the doubling of the number of peaks in the strong 
interaction regime is due do the fact that time reversal 
symmetry does not allow for one particle backscattering: 
one particle backscattering would lead to the expected 
Q peaks in pq(x,0) and ^(a;). The scenario can be in- 



FIG. 1. (Color online) Classical electron density po(a:,0) 
(panel (a)) and density-density correlation function ^(a;) 
(panel (b)), for Q=5 and Q/mL = 0.15. 


terpreted as a Wigner oscillation of quasi-particles with 
charge e/2: in the very same way strong interaction fa¬ 
vors Wigner oscillations of electrons with integer charge 
in usual spinful Luttinger liquids, here the interplay of 
spin momentum locking, two-particle backscattering, and 
electron-electron interactions, lead to the formation of 
such fractional oscillations. A physical insight in the 
meaning of this Wigner oscillation of fractional charges 
can be gained by refermionization: this technique allows 
to exactly solve the Hamiltonian in Eq.(l) for = 1/4, 
and for that value of the interaction parameter one has 
half charged fermions as low energy excitations. Since 
two-particle backscattering becomes relevant in the RG 
sense for Kl < 1/2, and we are considering Kl 
the quasiparticles of our system are also expected to be 
fermions with charge e/2. 

In order to refine the classical result, a theory for quan¬ 
tum fluctuations r](x) must be addressed. The appropri¬ 
ate framework is the collective coordinate method [50], 
and the most appropriate way of including fluctuations 
is to introduce them by the relation[49] 


/)(a:) = warn 
P 


2K(k)[^ + ^).k 


L 


2tt 


(9) 


which reduces to Eq.(4) for r]{x) = 0, up to the variable 
X. This behavior suggests the existence of a zero en¬ 
ergy mode rj(f{x) which is just a constant and stems for 
translational invariance. The expression for the electron 
density p{x) resulting from Eq.(9) can be conveniently 
written by using a standard series expansion for the Ja¬ 


cobi dn function as [51] 


p{x) = 


47r 


dxVix) 


OO 


n— — oo 


pin{4:'KQx I L+I3r}) 


cosh(nr) 


( 10 ) 


with T = 7rAr(\/l — k‘^)/K{k). The expression in Eq.(lO) 
closely resembles Haldane’s expansion for the electron 
density in the spinless Luttinger liquid [24], though with 
a striking difference: the harmonics of the density appear 
with wavevectors which are multiples of Akp = A’KQxjL^ 
instead of multiples of 2kp. This behavior is consistent 
with the oscillations in the density corresponding to the 
classical solution (Fig.l). Note that both in the limit 
m —)■ 0 (no cosine term) and in the limit Q/{LM) —>■ oo 
(significantly away from the Dirac point, that is where 
the effects of two-particle backscattering are expected to 
be negligible) only the term with n = 0 is present due to 
the presence of the damping factors cosh(nT). The usual 
form of the electron density of the hLL, i.e. the long wave 
density of usual one-channel Luttinger liquids, is hence 
recovered in such limits, as expected. 

The energetics of the fluctuations r]{x) is encoded, up to 
the quadratic order in r]{x), in the Lagrangian density[49] 

^riix) = [( 9 * 77)2 - [dxVp ( 11 ) 


with 


w{x) 


AK(k) AQK(k)x 

-dn --- ,k 

TT L 


( 12 ) 


We could exactly solve the equation of motion, de¬ 
riving from Eq.(ll), for the eigenmodes r]{x^t) = 

e'xp[—iujjt]r]j{x). In fact this equation becomes 

^jVzix) = w~'^{x)dx {w'^{x)dxVj{x )), (13) 

and it maps onto a Lame equation. The lowest en¬ 
ergy eigenfunction, which corresponds to w = 0 , is 
770 (x) = const., as required by the translation invari¬ 
ance of the system. When quantizing the theory this 
mode can be treated within the collective quantization 
approach, even though in the large Q limit it has been 
shown to be sufficient to integrate over 770 [49]. At low 
energy, still in the large Q limit, the theory for the low 
energy modes Pjix), with j ^ 0 , is essentially equivalent 
to a one channel Luttinger liquid [49], with Fermi veloc¬ 
ity uf = -s/l — k‘^K(k)/E{k) and Luttinger parameter 
K'^ = TTP^/{16y/l^K^{k)). 

By virtue of Eqs.(lO) and (11) in its low energy limit, 
we could show that the average electron density is, as 
expected, constant, and the zero temperature density- 
density correlation function ^(x) is, in the limit L —> 00 
and up to the slowest decaying oscillating term, 

^ L2 + (a;) (14) 
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with ^ 4 fcj^(a;) ~ cos ^ ^ j. The absence of the 

2kp component, which is the one that indicates the for¬ 
mation of a Wigner crystal, or at least Wigner oscilla¬ 
tions, in one channel Luttinger liquids, witnesses that 
the strongly interacting sector of the hLL is profoundly 
different from the usual strongly interacting electron gas, 
and that the hLL does not undergo a transition to a reg¬ 
ular one dimensional Wigner crystal. In fact two particle 
backscattering leads to a state characterized by charge 
fractionalization and the emergent quasiparticles (with 
fractional charge e/2) are hence 2Q due to charge con¬ 
servation. They interact among each other to form the 
fractional Wigner oscillations. 

In order to deepen the understanding of the effects of 
the original helical character of the theory on the strong 
interaction regime, it is natural to investigate how much 
spin momentum locking is affected when strong two- 
particle backscattering is included. A natural way to 
gain insight into this topic is to address the form of the 
spinor representing the Fermi operator. To do so one has 
to build the field 0{x), conjugated to rj{x). Leaving aside 
the discussion of the zero energy component, that is not 
relevant for the following, it is easy to show[49] that the 
field 6{x) = /^n(x')da;', with n(a;) = dC/d{dtr]) and 
with the quantization condition = iS{x—x'), 

can be evoked to properly define the Fermi operator for 
right/left moving electrons. By using the series expan¬ 
sion e““PK(/c)z] ^ ^ 

TYgTi'n-+'^/'^)/(^kK{k) smh{2n + 1 )t), and the definition of 
the chiral fields (j)R/L{x) = j/)(a;)=F we obtain the 

expression for the Fermi fields in terms of r]{x) and 0{x) 

OO 

(-a)g±z(4^Qa:/L+/3r,)(n-H/4)^ 

n——oo 

(15) 

where Cn '^ is implicitly given by the relation = 
■ The form of the Fermi spinor 4'(x) 
is hence 4'(a::) = {tp+{x),ip-{x))'^ = ('(/^(a:),V' l(x))^, 
where ip±{x) is the spin up/down component of the field 
operator. The usual form of the hLL can be recovered 
by setting Cn^ oc as can be easily shown to be 

the case in the limits m —>■ 0 and Q/{mL) —>■ oo. For 
finite m and Q/(mL) the breaking of spin momentum 
locking is evident: the terms with n > 0 (n < 0 ) 
represent left (right) moving components in the spin 
up (down) part of the spinor. However these terms are 
suppressed by both the coefficients Cn'^ and the scaling 
they acquire in the correlation functions. Since it is a 
direct consequence of the chirality of the Fermi spinor, 
the strong anisotropy in the spin response function, 
hallmark of the hLL, is hence to some extent also 
present in its strong interaction sector in the presence 
of two-particle backscattering. This is in sharp contrast 
with the behavior of the one dimensional Wigner 


crystal, which is an almost classical state in which 
spin dynamics plays an unimportant role, and is often 
integrated out [42]. To support our claim we address 
the ground state average Sij{x) of the spin-spin func¬ 
tions Sij{x) given by Sij{x) = 'k^'(a;)(T*'I'(a:)'I'l(0)(T^4'(0). 

In the hLL without two-particle backscattering one 
finds that sii(a:) = S 22 (x) and 512 ( 2 :) = —S 2 i(a:) are 
purely oscillating functions with wavevector 2kp, while 
533 ( 2 :) is a non oscillating function and the remaining cor¬ 
relations vanish. In our case, 533 ( 2 :) is still a non oscillat¬ 
ing function, so that the probability of finding two elec¬ 
trons at distance x ^ a from each other is independent of 
the spin projection and follows the density-density corre¬ 
lation function. On the other hand the oscillations in 
511(2:), 522(2:) and 512(2:) have a wavelength which is 
twice the wavelength of the oscillations of the electron 
density and hence they witness the onset of a spin helix: 
this behavior is in accordance with what happens in the 
weakly interacting hLL. However the spin oscillations in 
the strongly interacting edge with two particle backscat¬ 
tering are built on oscillating density-density correlations 
characterized by half of the wavelength, suggesting an in¬ 
tercalation of electrons with spin up and spin down, once 
a projection axis perpendicular to the 2 ; axis is chosen. 
An intuitive picture is drawn in Fig. 2. 



FIG. 2. (Color online) Scheme of the strongly interacting 
ground state: the fractional quasiparticles tend to form a 
Wigner crystal (the density-density correlation function is 
shown). Moreover the spin density is schematically depicted 
by using the arrows. 


In conclusion we have shown that two-particle 
backscattering, in combination with strong electron elec¬ 
tron interaction, leads to the formation of a peculiar cor¬ 
related structure, a Wigner crystal of quasi-particles with 
charge e/2. This state still bares the signatures of the un¬ 
derlying helicity of the system, both in the form of the 
Fermi operator and in the anisotropy of the spin-spin cor¬ 
relation functions. Moreover such a state has a complex 
spin structure, resembling a spin helix, that suggests an 
intercalation of particles with opposite spin projection 
once a projection axis perpendicular to the z axis is cho¬ 
sen. 
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